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Abstract Relativistic shocks provide an efficient method 
for high-energy particle acceleration in many astrophys- 
ical sources. Multiple shock systems are even more ef- 
fective and of importance, for example, in the internal 
shock model of gamma-ray bursts. We investigate the 
reacceleration of pre-existing energetic particles at such 
relativistic internal shocks by the first order Fermi pro- 
cess of pitch angle scattering. We use a well established 
eigenfunction method to calculate the resulting spectra 
for infinitely thin shocks. Implications for GRBs and rel- 
ativistic jets are discussed. 
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Relativistic Shocks 



r— j • 1 Introduction 



Observations of the relativistic jets associated with ac- 
tive galactic nuclei (AGN) show strong non-thermal emis- 
sion at different sites within the jet structure and where 
the jet meets the intergalactic medium. While we can- 
not observe an outflow in gamma ray bursts (GRBs) a 
similar, but much more relativistic, jet structure is as- 
sumed with the prompt emission coming from the inter- 
nal shocks while the external shock is responsible for the 
afterglow. 

By analogy with the non-thermal emission from the 
non-relativistic shocks in supernova remnants, a first or- 
der Fermi process has been proposed by various authors 
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as the mechanism responsible for the production of non- 
thermal particles at relativistic shock fronts. The semi- 
analytic eigenfunction method first introduced in the late 
80's [H[5] allowed for calculation of the spectral index and 
anisotropy at shock front moving up to Lorentz factor 5 
in the upstream rest frame. An improvement in Kirk et 
al. [4] allowed a rapid computation of the spectral index 
of arbitrary shock speeds with the help of an analytic ap- 
proximation. The spectral index of particles accelerated 
at such shock fronts has been shown to tend to 4.22 for 
increasing Lorentz factors, in the hydrodynamic limit, 
where the compression ratio tends to 3. It has also been 
shown that the particle distribution is anisotropic at the 
shock front even for mildly relativistic flows. 

Ignoring the injection mechanism, the jet produces a 
power law spectrum of electrons at a hydrodynamically 
dominated external shock via first order Fermi accelera- 
tion. This natural spectrum of the external shock, s, is 
calculated using the method of Kirk et al., of which we 
have given a brief summary. We then let these particles 
escape downstream where they are accelerated again by 
an internal shock, which has natural index q > s. Adapt- 
ing the work of Kirk et al., we are able to calculate the 
number of particles at the internal shock front compared 
to the number far upstream at the external shock. We 
can also measure the anisotropy at the shock and the 
number of particles far downstream of the internal shock. 

In Section 2 we summarise the semi-analytic eigen- 
function method for determining the spectral index of a 
power law distribution produced at a shock and derive 
the distribution function far downstream of the internal 
shock. In Section 3 we describe the conservation of par- 
ticle flux and use it to motivate a consistency check on 
our results, which are presented in Section 5. 



2 An Eigenvalue Approach to The Transport 
Equation 

We restrict ourselves to the acceleration of particles with 
momenta well above the injection momentum but below 
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any cut-off momentum where radiative fosses change the 
shape of the spectrum, fn the shock rest frame the parti- 
cie transport equation, which holds separately upstream 
and downstream, can be written as 

where u is the fluid velocity, r = (1 — it 2 ) -1 / 2 , D^^ is 
the pitch-angle diffusion coefficient and / is the parti- 
cle phase-space distribution. We have to ensure that the 
distribution matches at the shock. 

/-(p_,M-,0) = /+(p + ,/x+,0) (2) 

where the plus (minus) sign denotes quantities upstream 
(downstream) of the shock. 
We can expand / as 



f(p,IJ>, z)= &iO)Qi(M)exp ( 



i— — oo 



where (Qi(fj,), Ai) are an eigenfunction, eigenvalue pair 
satisfying 

^ D w^ = Mu + »)Qi (4) 

and the eigenvalues are ordered such that 

A^i-i < A_i < A — < Ai < A i+ i Vi > 

To find a bounded solution with no particles far up- 
stream we set bi(p) = ciip~ s and the upstream solution 
becomes 



/ = z2 a * P-'Qi (^-) ex P 



i>0 



a: 



r 



(5) 



with normalisation a 1 =1, while the downstream can 
be written as 



/+ = 5>+p+ s Q+0 U+ )e X p 



i<0 



A 



r 



(G) 



then the matching condition can be reduced to 

Yj^P-Q-in-) = Y J 4v+Qt^+) (7) 



i>0 



i<0 



Using the orthogonality of the eigenfunctions we can 
then solve for s and the a^s as in Kirk et al. In the 
limit of non-relativistic this gives the usual s = 4 result 
with negligible anisotropy, while for highly relativistic 
shocks we reach s = 4.22 with considerable anisotropy, 
approximately described in the shock rest frame as 

f(p s h, Msh, 0) oc p^(l - /i sh u-)~ s exp ( ( 1 + ^ sh . 

Now suppose far upstream we have a known particle dis- 
tribution g{jpJ). Then our matching condition for bounded 
solutions becomes 

g(p-)+YbT(p_)QT (,,_)= Ybf(p + )Qt(v+) (8) 



i>0 



If our far upstream distribution is a power law g(p-) = 
p_ q we can find solutions with bf(p±) ~ afp± q . 

Using the Lorentz transformation p + = _T rc ip_(l + 
u rc i/i_) where w ro i = (m_ — u + )/(l — u_w + ) the matching 
condition reduces to 

(1 + UrolM-)" + + "rclA*-) 9 Qr(/ i -) 

i>0 

= Y a t r JQi(»+) (9) 

i<0 

Then multiply by (u++/j+)<3+(/i+), j > 1, and integrate 
over to get 

^aZW id = -W OJ /Qo (10) 

i>0 

where 



(3) 



(l+u IcU ^) q (u + +^ + )Q+{^+)Q l ( A i_)d//+(11) 



so in matrix form a = W _1 b where bj — —Woj/Qq. 

We now have the solution at the shock in the up- 
stream rest frame. The isotropic distribution far down- 
stream is thus 

1 f 1 

MP+) = o — / / (P- 1 M-)S = 0)(tt++jii+) ef/i+(12) 
2u+ J_i 



3 Conservation of Particle Flux 

In Heavens & Drury [5] a consistency check was intro- 
duced to ensure conservation of particle flux in phase 
space. It can be explained as follows: consider all parti- 
cles (upstream and downstream) with momentum up to 
p+ as measured in the downstream rest frame. In the 
upstream rest frame this is all particles up to p- = 
r re \p+(l + u re \). The amount of particles in this range 
can only change, in the steady state, by, injection at mo- 
menta lower than p + , flux from far upstream, flux lost 
far downstream and by particles accelerated at the shock 
by crossing from downstream to up and returning with 
momentum p* + = + it rc i)/(l — u re i/i+)- This can be 
written as 

rl rP- 

<P +r_ / / (u-+fX-)f(p'_,ix-,-oo)2Trp'ldp'_dfJ,- 



= r + 



-1 Jo 



1 Jo 

1 r*% 



(u + + fi + )f(p' + ,fj, + , oo)2Trp+dp' + dfj, A 



i<0 



+ r+l I (u + +^ + )f(p' + ,f, + ,0)2np^dp' + dfx + (13) 
J -i Jp+ 

where <P is the integrated injection flux up to p+. Given 
/(p_,/i_, — oo) = pZ q , we differentiate with respect to 
p + and arrive at 

2r_u_p- q (r Icl (i + u IC ,))- q + 3 

-r+ ( u + + M+)/(P+>M+,°o)<^i + = 
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(1 - M rc l/i+) 3 

/(p+,A*-H )W+ 



(14) 

Assuming the solution isotropises far downstream, the 
obvious solution is f(p + , fi + , 0) = p+ q g{^+), f(p+, (J.+ , °o) 
goP+ q where go = J (u + + n + )g(n + )diJ, + /2u + and equa- 
tion (fl"4f reduces to 

2 J r_ M _(r rc i(i + Mrc i)) 3 -9 
(i + Urcl ) 3 -" 



= r+ / (u + + fi + ) g(fJ,+) T rg— d/i+ (15) 

J-i V (i - ur C in + y v 

When g = 3we see that go = R, the proper compression 
ratio. For q ^ 3 we use 

., r + f\{u + + (i + )g(ji + )(l - u rcl [i + ) c i- 3 dti + 
JVc = 5—; (16) 

as a measure of the consistency of our results, where 
A/" c = 1 implies an exact solution . 



4 Adiabatic Gains 

In the absence of a diffusion process capable of scatter- 
ing downstream particles back upstream, particles will 
undergo an energy change due to the compression of the 
plasma. Since the position element of phase space, d 3 x 
is compressed by R, the momentum element, d 3 p, is ex- 
panded by R. Hence the downstream momentum can be 
related to the upstream momentum by p + = i? _1 / 3 p_. 

-oo) = p~ q then the adiabatic 
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Fig. 1 The effect of magnetisation on the amplification pro- 
cess with a p~ 4 spectrum far upstream. As the magnetisation 
increases the compression ratio decrease, as does the ampli- 
fication. On the x-axis we have the 4-speed of the shock as 
measured in the upstream rest frame. 



First we will define a magnetisation value 



1 



(18) 



Thus if we have /(p-,- 
gain is 

! ( {Z - °°\ = R^ (17) 

This is often referred to as "shock-drift" acceleration but 
it is in fact only the simplest form of it (for a more de- 
tailed view of shock-drift see Begelman & Kirk pQ). We 
use this as a benchmark as it is the least amplification a 
power law can under go at a shock. 



5 Results 

The results presented in this poster consider the back- 
ground plasma to have perpendicular magnetised shocks 
with negligible upstream hydrodynamic pressure, although 
it should be noted that they are also relevent for hydro- 
dynamical shocks where the upstream pressure is impor- 
tant, resulting in relatively weak shocks. This would be 
the case if the internal shock's upstream was the down- 
stream of an external shock. The equation of state is 
taken to be Juttner-Synge with the downstream quanti- 
ties solved from their upstream counterparts as in Kirk 

6 Duffy [3]. 



where va is the Alfven speed. 

Figure Q] illustrates the additional energy gain ob- 
tained from diffusive shock acceleration over shock drift 
acceleration for a population of cosmic rays with a power 
law distribution p~ A advected into various relativistic 
shocks. Care must be taken when interpreting the re- 
sults in the hydrodynamic case for low ru as in this 
case s —> q = 4 and the contribution of injected parti- 
cles will be very important. However for large ru the 
advected cosmic ray population with p~ 4 will dominate 
over the p -4 22 population originating from injection at 
the shock. It should be noted that the gain is more than 
one order of magnitude greater than that obtained via 
shock drift acceleration. So when ever the downstream 
plasma expands and the shock disappears the cosmic ray 
population will have been genuinely accelerated as the ef- 
fects of adiabatic losses will not cancel the energy gain of 
the diffusive process as it would have the gain by shock 
drift. 

Shocks with high a and a magnetic field perpendic- 
ular to the velocity have a reduced compression ratio, 
leading to an increase in their natural spectral index. 
However if there is a power law with a harder index 
upstream then we will still have amplification of it by 
the magnetised shock, and while it may be smaller than 
the amplification by its hydrodynamic equivalent, it is 
several orders greater than the non-relativistic predic- 
tion. Again the energy gain is greater than that achieved 
by shock drift acceleration, although we can see from 
the example of a = 0.1 that the energy gain difference 
is small for mildly relativistic shocks. However as the 
shock becomes more and more relativistic the difference 
between the energy gain from the two different processes 
increases. In fact, by looking at figure[2]it seems that the 
energy gain from the diffusive process goes like r^T 2 for 
ultrarelativistic shocks, while the energy gain for shock 
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Fig. 2 We measured the amplifications in the downstream 
rest frame for three power law distributions far upstream, 
p~ 4 , p~ 41 and p -4 ' 2 . On the x-axis we have the 4-speed of 
the shock as measured in the upstream rest frame. It seems 
that the energy gain from the diffusive process goes like r^~ 2 
for ultrarelativistic shocks. 



drift is R^ 3 
we have 



r 



q/3 



Thus for ultra-relativistic shocks 



energy gain DSA 
energy gain shock drift 



oc r 



2(g-3)/3 



(19) 



This is the same scaling as noted in Begelman & Kirk 
when they took a more detailed look at "shock-drift" 
acceleration in the ultrarelativistic limit. 



5.1 Accuracy and Consistency of the Method 

While the method has been stated in a purely analytic 
manner, finding our cigcnfunctions requires numerical in- 
tegration, see Kirk et al. for details. As the left hand 
side of equation ([3]) changes sign at fi — —u accurately 
sampling the functions between fj, = — 1 and /i = — u is 
crucial, especially for highly relativistic flows. We chose 
a simple non-adaptive step size in [i, Afi, and checked 
for the convergence of our solutions as we reduced Afi. 
The results shown in figure [3] are for a power law with 
index q = 3.95 advected into an unmagnetised strong 
shock satisfying the Juttner-Synge equation of state. 

Figure [3] shows how the convergence of the computed 
far downstream distribution with decreasing A\i. As ex- 
pected the larger the Lorentz factor of the upstream 
medium in the shock rest frame the smaller we require 
A/i to ensure an accurate result. Also the measure of 
consistency, Af c , equals 1 for small enough Z\/i. 



6 Application to Jets with Internal Shocks 

If the shock is propagating into the ISM or IGM then we 
would expect a upstream of the shock to be small. If the 
shock is relativistic then it is expected that it is also a 
perpendicular shock. For perpendicular shocks it can be 
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Fig. 3 Convergence of the far downstream amplification, 
where upstream we have p~ 3 ' 95 , as A/i decreases. In our 
range of interest, F < 100 the results converge completely 
by An = 4 x 10~ 7 . 



shown by a simple manipulation of the results presented 
in Kirk & Duffy that 



+ 



1 



(20) 



1t+ \1 + <T_(1 — U-/U+), 

where the velocities are measured in the shock rest frame. 
For small <r_ at ultrarelativistic shocks (u_ = 1) v+ « 
1/3 and 

a+ w 3cr_ (21) 

Thus the downstream of such a shock can have a high 
a = <r5_, unlike parallel shocks which reduce a. Due to 
the Lorentz invariance of a the upstream magnetisation 
value of this shock a % _ will equal er^ . This high magneti- 
sation value reduces the compression ratio of the second 
shock allowing it to be a hydrodynamically weak, rela- 
tivistic shock. Such shocks are not capable of producing 
hard power laws by themselves. However, assuming the 
external shock produces a hard power law, they can act 
as effective reacceleration sites throughout the jet, using 
the same pitch-angle diffusion mechanism responsible for 
production of hard power laws at strong shocks. 
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